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ABSTRACT 

We report Submillimeter Array observations of SO2 emission in the circum- 
stellar envelope of the red supergiant VY CMa, with an angular resolution of 
~ 1". SO2 emission appears in three distinct outflow regions surrounding the 
central continuum peak emission that is spatially unresolved. No bipolar struc- 
ture is noted in the sources. A fourth source of SO2 is identified as a spherical 
wind centered at the systemic velocity. We estimate the SO2 column density 
and rotational temperature assuming local thermal equilibrium (LTE) as well 
as perform non-LTE radiative transfer analysis using RAD EX. Column densities 
of SO2 are found to be ~ 10^^ cm~^ in the outflows and in the spherical wind. 
Comparison with existing maps of the two parent species OH and SO shows the 
SO2 distribution to be consistent with that of OH. The abundance ratio fsOi/ fso 
is greater than unity for all radii greater than at least 3 x 10^^ cm. SO2 is dis- 
tributed in fragmented clumps compared to SO, PN, and SiS molecules. These 
observations lend support to specific models of circumstellar chemistry that pre- 
dict fso2/fso > 1 and may suggest the role of localized effects such as shocks in 
the production of SO2 in the circumstellar envelope. 

Subject headings: stars: individual (VY CMa) — stars: late-type — circumstellar 
matter — submillimeter — radio lines: stars 
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Introduction 



Massive stars (M^, > SM©), enter a red supergiant (RSG) p hase during which th e 



star experiences mass-loss at rates of M ~ 10 ^ — 10 ^MqUt ^ (Ivan Loon et al 



20051) 



he time variation of this mass-loss rate is not well-constrained by theoretical studies 



( lYoon and Cantiello 



20101 ) ■ As a result, the total amount of mass lost over the course of 



the RSG phase remains uncertain for a given initial mass ( ISmith et al. 



20091 ) ■ Observations 



of ma ss-loss events have shown them to be sporadic and spatially anisotropic (ide Wit et al. 



20081 ) ■ 



VY Canis Majoris (VY CMa) is an oxygen-rich red supergiant with an estim ated mass of 
25 M(r, and a mass -loss rate estimated to be M ~ 2 — 4 x 10~^MQ?/r~^ flDanchi et al. 



1994 



Smith et al. 



20091 ). Optical images of this source show multiple discrete and 



asymmetric mass-loss events, ranging in a ge from 1700 to 157 yea rs ago, that are distinct 



from the general flow of diffuse material (IHumphreys et al. 



20071 ). Detailed studies of 



mm/sub- mm molecular spe ctra of VY CMa haye been carried out, revealing the chemica l 



complexity in the envelope (IZiurys et al 



2007 



Royer et al 



2010 



Tenenbaum et al 



Spatial structures in the visible and IR bands have also been obtained (iSmith et al. 



High angular resolution Submillimeter Array (SMA£ 



maps of the spatial distribution of CO and SO (IMuller et al. 



201G) 



200 ih . 



observations of VY CMa produced 



20071). 



Millimeter wavelength observations have long show n evidence of SOo in th e circumstellar 



envelopes (CSE) of oxygen-rich red supergi ants (e.g., 



emission in VY CMa was first detected by 



Guilloteau et al. 



Sahai fc Wannierl (jl992l ). Later infrared 



19861) and SO2 



observations have also shown the production of SO2 within several radii of 0-rich AGB 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded 
by the Smithsonian Institution and the Academia Sinica. 
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stars (in the "inner wind") lYamamura et al.l (Il999l ). From a theoretical perspective, the 
presence of SO2 in circumst ellar envelopes has bee i i predicted by isotrop i c, non-equihbriurn . 



models of stellar chemistry (IScalo &: Slavsky 



1980 



Neiad fc Millar 



1988 



Willacv fc Millar 



19971 ). This class of models assume an isotropic geometry of the CSE with specified values 
for the mass loss rate and expansion velocity. The resulting model CSEs consist of an 
inner region where assumed parent species in the outflow are broken down, a intermediate 
region where chemical reactions lead to high abundances of daughter species, and an outer 
region where the daughter species are destroyed by photodissociation. SO2 in these models 
has been assigned as a daughter species, although it has since been s hown to exi s t as a 



2010h . ICherchnefB (|2006|) 



possible parent species created in the inner wind (iDecin et al. 
studied the role of shocks in SO2 formation in the inner wind within a few stellar radii of 
the photosphere. 



Jackson & Neuven 


1988; 


Yamamura et al. 


1999; 


Decin et al. 



20101 ) have repeatedly shown SO2 abundance to be higher both in the ex tended CSE and 



the inner wind than ex pected from CSE and inner wind chemistry models (iWillacy fc Millar 



1997 



Cherchnefl 



passage of shocks (IHartquist et al. 



20061). Local abundanc e of SO2 may be enhanced in the ISM by the 



19801 ). and shock chemist ry has been proposed to 



explain the high observed SO2 relative abundances in C SEs (IJackson fc 



Alternativ ely, low model a 



to that of 



j undan ces of SO2 presented in 



Scalo &: Slavsky 



Nguyen 



Neiad fc Millarl f ll988f ) compared 



19881). 



(Il980l ) may be due to the latter's lower assumed value for the 



photo dissociation of SO;^, wh ich is the main mode of destruction of this molecule at large 



radii (iSahai fc Wannier 



19921 ) 



Interferometric mapping of SO2 distribution in the CSE of an evolved star may 
contribute to the improved understanding of sulfur chemistry by directly constraining 
the relative abundances of sulfur-bearing species as a function of radius. Furthermore, 
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associations of SO2 enhancement with discrete outflow features may favor the idea 
of locahzed production a nd provide evidence for non-isotropic processes in CSEs 



( jjackson fc Nguyen 



19881). In this paper we present results of high spatial resolution maps 
obtained from SMA observations of SO2 around VY CMa. We describe four spatially 
discrete sources of SO2 emission and use multiple transitions of SO2 and derive rotational 
temperatures and column densities (assuming local thermodynamical equilibrium). 
Furthermore, we perforin non-LTE radiative transfer analysis using the RADEX package 



( Ivan der Tak et al. 



20071 ) to derive local kinetic temperatures, H2 densities, and SO2 



abundances of the SO2 emitting regions. We then compare our maps to published visible, 
IR, and submillimeter observations to evaluate their consistency with the above cited 
models of circumstellar chemistry. 



2. Observations 

VY CMa was observed with the SMA on 2009 February 18, with the array in extended 
configuration, offering baselines from 44.2 to 225.9 m. The projected baseline lengths 
were from 14 to 188 m. The frequency coverage was 234.36 to 236.34 GHz in the lower 
sideband and 244.36 to 246.34 GHz in the upper sideband. The phase center was at 
a(2000) = 07'^22"^58^27, 5(2000) = -25°46'03.4". The quasars 0730-116 and 0538-440 were 
observed every 20 minutes for gain calibration, and the spectral band-pass was calibrated 
using quasar 3c273. Flux calibration was done using recent SMA measurements of 0730-116 
(2.75 Jy at 1 mm) and 0538-440 (3.96 Jy at 1 mm). Nominal flux calibration accuracy is 
15 to 20%, depending on the phase stability. In our observations, the uncertainty appears 
to be better than 10%, based on the good agreement with the ARO spectra (see section 
3.2 and Figure 3). The on-source integration time on VY CMa was 5.66 hours. T^^^ (SSB) 
varied approximately from 200 to 400 K during the track with an atmospheric zenith optical 
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depth of ~0.1 at the standard reporting frequency of 225 GHz, measured at the nearby 
Caltech Submilhmeter Observatory. The conversion factor between Kelvin and Jansky for 
our observations is ~ 14 K/Jy. 

The visibihty d ata were cahbrated using the M IR package in IDL and imaged with the 



MIRIAD softwarcl flSault. Teuben fc Wright 



1995|). The field of view (FWHM primary 



beam) varies from 53" to 56" whereas the largest angular extent of the source is expected 
to be about 5". The synthesized beam size, representing the obtained spatial resolution, is 
l."48 X 1.03". The adopted weighting mode for imaging is "natural" weighting (specified 
in the Miriad task INVERT). The final rms noise level is ~60 mJy/beam per channel in 
the spectral line images and 3.9 mJy/beam in the continuum map. Line emission was 
subtracted b yvisually examining the spectra in visibility amplitudes, and using the line-free 
channels specification to the Miriad task UVLIN. 



3. Results 

3.1. Continuum emission 

We detected continuum emission in both the upper and lower sidebands, which is well 
described by a point source centered at RA=07''22"'58.336', Dec=-25°46'03.063". A 2D 
Gaussian fit in the image plane yields an integrated fiux of 335.0 ± 66 mJy at 235.4 GHz 
and 359.3 ± 72 mJy at 245.4 GHz (maximum errors are of 20%, mainly resulting from 
the uncertainty in fiux calibration). Previous SMA observations of VY CMa at 215 and 

uxes o f 270 ± 40 mJy and 288 ± 25 mJy respectively 



( IShinnaga et al. 



2004 : 



continuum f 


uxes c 


MuUer et al. 


2007) 



20071 ). These four values are consistent under the 



Rayleigh- Jeans approximation with black-body radiation with a brightness temperature 



^http:/ /www. cfa.harvard.edu/~cqi/mircook. html 
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higher than 12 K. 



3.2. Line emission 



In the two sidebands of ~4 GHz bandwidth, we detected a total of 14 hnes, which show 
emissions from CS, H2O, PN, SiS and SO2. In Figured] we present the spectra and maps 
of the SiS J=13-12 and the PN J=5-4 lines. These lines are representative of a relatively 
compact emission that is centered on the peak of the continuum emission, hence they allow 
us to estimate the systemic velocity of the star. The SiS s pectrum shows the distin ctive 



triple peak morphology as described by 



Ziurys et al. 



(120071 ) 



Tenenbaum et al. 



mm that is 



interpreted as the signature of a slowly expanding shell and a pair of faster outflows nearly 
collimated with the line of sight. We use the fitted velocity of the central peak of these lines 



to find the systemic velocity of the star vlsr = 19 .5 km s ^. This va . 



other results from mm and sub-mm observations (IZiurys et al 



2007 



ue is consistent with 



MuUer et al 



2003) 



and is adopted as the star rest-frame velocity in the following sections. 

For the remainder of this work we focus only on the SO2 lines, which are summarized 
in Table [H The excitation energies of the four identified SO2 line s vary from 19 to 130 K . 



Tenenbaum et al. 



torn 



All four of these lines were detected in the spectral line survey of 
Our interferometric observations reveal the spatial distribution of the SO2 emission for the 
first time. The spectra and integrated maps of all four lines are shown in Figure [2J Figure 
|3] compares the spectra of two SO2 lines with single dish data from the Submillimeter 
Telescope of the Arizona Radio Observatory, while channel maps of all four lines are 
presented in Figures HI El [6] and [71 

Maps of SO2 lines maps show well resolved spatially extended emission over an area of 
3" X 5", in which we distinguish four distinct components, hereafter designated as sources 
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A, B, C, and D. Figure |8] shows all four sources in a position-velocity space diagram of 
the 245563.4 MHz line. Source positions, orientations, and intensities were fitted to 2D 
Gaussians using the MIRIAD routine IMFIT. Source A is a very strong blue-shifted source 
with visr between —18 and +10 km s~^. It is spatially offset to the East of the stellar 
position by ~ 0.7". Of similar intensity to source A is source B, an elongated source offset 
by ~ 0.8" to the West. Source B has a broader spectral profile than source A and is heavily 
red-shifted with visr between +26 and +66 km s~^. Source C is a strong source with a 
similar radial velocity as source B. In contrast to sources A and B, however, it is heavily 
offset to the West by - 2.5". 

Source D is significantly weaker than sources A, B, and C. In velocity space, it is 
centered on the star's frame of motion (12 km s^^ < Vlsr< 26 km s^^) with a distinctly 
narrower spectral profile and similar spatial dimensions as the other sources. It is also 
spatially centered at the star's location. 

Assuming local thermal equilibrium, we can use the integrated intensities measured in 
the detected SO2 lines in each of the four identified SO2 sources to constrain t heir respective 



column density and rotational temperature Trot, using the expression from 



Friedel et al. 



fl2005h : 



beam 



In 



Eu 



■IGir^Beaebiy^Sjj,'^' '"'Z' Trot 
where heam is the integrated flux in Jy km s~^ beam^^, 6a6b is the beam size, and B 



from 


Snvder et al. 


(2005) 


(Miiller et al. 


2005 


), S is 



(I2OO5I ): Q source /{Q beam + Q source)- Z is the rotatioual partition function 



20051 ). S is the line strength, /i the dipole moment (1.63 Debye for SO2), Eu 



the line's uppe r rotational state energy in K taken from the JPL spectroscopic database 



(Pickett 



199ll ). and 9 the solid angle. 
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By performing linear fits of the rotational temperature diagram shown in Figure |9l 
we find that the three faster outflows (sources A, B, and C) show variations in rotational 
temperature. Sources A, B, and C have rotational temperature of Gl^f^ (one a error), 
110^29 69127 respectively. The column densities (N) of all three sources are 

similar (A: 9.2tli x 10^^ cm'^; B: l.Ot'^oi x 10^*^ cm^^; C: l.ltol x 10^^ cm^^). The 
rotational temperature of the close-in source D is much higher at 240^52 K, while its column 
density is likely lower than that of the others at 7.6 ± 1.0 x 10^^ cm~^. 

The derived parameters for source D are uncertain since a significant portion of flux in 
its velocity range was missed by our interferometric observations due to the lack of short 
spacing data (our shortest projected baseline 14 m). Comparison with single dish data 
from the Arizona Radio Observatory shows that features as large as 4" are fully recovered 
(Figure [3]). Features much larger than 4" may be subject to an underestimation of its flux, 
although 4" is an highly conservative estimate. To address this deficiency, we repeated the 



rotational temperature diagram analysis for source D using a single dis 



1 observation made 



by th e Arizona Radio Observatory 10 meter Submillimeter Telescope (ITenenbaum et al. 



2010l ). Assuming a source size of 4" by 4" we derive a lower rotational temperature of 



97l3^ and a column density of 1.9l^;g x 10^^ cm^^. Because this column density represents 
an average over our assumed source size, it cannot be directly compared to our derived 
column densities, which represent the values at the center of each source. This averaged 
column density also represents an upper limit due to our assumption of the smallest possible 
angular size. 

The validity of the LTE approximation can be assessed by checking the linearity of the 
data points in the rotational temperature diagram. Deviation from a linear trend suggests 
non-LTE conditions in the source or line misidentification. In the case of source D, deviation 
from linearity is very small compared to the error. We are therefore confident that LTE 
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is a valid assumption for this source. For the outflow sources A through C, a systematic 
concave up shape is noted in the rotation temperature diagram, suggesting departure from 
LTE conditions. 

These findings justify the need to consider a non-LTE radiative transfer model to 
i nterpret the measured f luxes. We perform this analysis using the RADEX package 



( Ivan der Tak et al. 



20071 ). which calculates expected line intensities for a given column 
density of SO2, kinetic temperature (Tkin), and volume density of H2. For each source, we 
assume spatially homogeneous kinetic temperature and volume density. The SO2 column 
density is fixed for each source and its value is drawn from the LTE analysis, which gives 
10^^ cm~^ for sources A, B, and C and 10^^ cm~^ for source D. Varying column density one 
order of magnitude in each direction does not significantly affect the results. We perform 
calculations for a wide range of values for the kinetic temperature and the H2 density. We 
then evaluated the resulting ratios between the flux of each line and that of the 235151.7 
MHz transition for the outflow sources A, B, and C. Due to the lack of emission from the 
235151.7 MHz line in the vicinity of source D, the 236216.7 MHz line was used instead to 
calculate line ratios. 

Fits to our data are presented as the statistic p-value for each assumed value of 
Tkin, and H2 density (Figure [101) and best fit results are tabulated in Table [3l We see that 
the non-LTE RADEX results for Tkin show general agreement with Trot derived with LTE 
assumption. A reliable estimate of H2 density is elusive. Because the flux of each line 
becomes more similar at greater values of H2 density, we are able only to constrain a lower 
bound on this value. 

The constraints on the temperature and H2 density in source D are much weaker. 
However, the best fit values of Tkin span the 130 K to 330 K range, which brackets the value 
of Trot = 240 K obtained from the LTE analysis. This agreement between Trot and Tkin is 
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consistent with LTE conditions in source D. Our RADEX analysis also allows us to evaluate 
our assumption of optically thin SO2 lines necessary for the LTE analysis performed above. 
Although for lower H2 densities (< 10'' cm~^), the 235.151 GHz line comes close to being 
optically thick, for the inferred H2 densities, all four SO2 have optical thickness below 10^^. 

Assuming that the extent of the sources in the line of sight direction is similar to that 
in the plane of the sky, and that the density is homogenous, we estimate the fractional 
abundance of SO2, /502 i^i sach source. These are tabulated in Table [31 



4. Discussion 



4.1. Source properties 



We begin by discussing the positions of the four identified SO2 sources. Source C, offset 



2.5" to the West, is clearly an isolated 
existence was suspected in IZiurys et al. 



3 ody w ith no antipodal companion to the East. Its 



( 120071 ). but was not treated (IS db distinct source. 



On the other hand sources A and B are found to be offset in opposite directions relative to 
the star position with similar blue and redshift velocities relative to the stellar frame. SMA 
observations of the CO and SO analogs of sources A and B were interp reted to be antipodal 



companions o riente d 15° from the line of sight by 



MuUer et al. 



Ziurys et al 



fcoOTn . On the other hand. 



(120091 ) found that single dish observations of CO and other molecules are best 



explained by a blueshifted and a redshifted source at 20° and 45° from the line of sigh t 



Similarly, visible and IR HST observations ( iHumphreys et al 



2007 



Smith et al 



2OO1I) 



have found no evidence of antipodal structure around VY CMa. A careful inspection of 
sources A and B in our observation shows that the faster (in radial velocity relative to the 
star's reference frame) sections of both bodies are offset towards the southwest, while the 
slower sections are offset towards the northeast. This observation argues against a bipolar 
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geometry, in which antipodal subsections of the two outflows should have similar velocities 
We therefore treat sources A, B, and C as three distinct outflows probably unrelated to 
the symmetry axis of the star. These sources do not seem to correspond to any visible 
or IR features. Our maps show that all SO2 outflow sources are too close to the star to 



be ide ntified as the "curved nebulous tail" or the numbered arcs presented in 



Smith et al 



( I2OOII ). which are located ~ 3.5" from the star. 



We i dentify sources A and B with the blu e and redshifted ou tflows of the previous 



authors fiZiurys et al. 



2007 



MuUer et al 



2OO7I 



Ziurys et al 



20091 ). The P-V space 



mo rphology of these so urces match closely with outflows of both CO and SO lines mapped 



by 



MuUer et al.l ( 120071 ) (Figure |8]). These bodies are hypothesized to have originated in 



an episode o 



( ISmith et al 



anom alously high mass loss at uncorrelated locations on the stellar surface 
20011). Ass uming that th e ages of the outflows are similar and adopting the 



500 year age found by 



MuUer et al. 



(120071 ). we can attempt to derive their locations in 
three-dimensions. We find that the deprojected radii of sources A and B are 4.2 x 10^^ and 
4.8 X 10^^ cm and that both are situated at 22° from the line of sight. Their deprojected 
star- frame velocities of 28 and 30 km s~^ fall within the range of measured outflow velocities 



from the multi-epoch observations of 



Humphreys et al. 



(120071 ). If we assume the same age 
as for sources A and B, then source C is found at a similar radius from the star (6.9 x 10^^ 
cm), but it is much faster at 44 km s~^ and is situated 54° from the line of sight. SO2 



abundance at large radii is expected to be co ntrolled by a balance 



Detween the rates of 



1980l ). Outflows with faster 



production, expansion, and photodissociation ( iScalo fc SlavskyI 
expansion velocity are expected to maintain high SO2 abundance out to greater radii, as in 
the case of our source C. 

Source D is elongated and centered at the stellar position, with a resolved minor and 
major radii of 1" and 1.6", corresponding to 2.2 x 10^^ and 3.5 x 10^^ cm. The minor 
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axis radius co rresponds to between 180 and 540 stella r radii, depending on the adopted 



stellar radius ( iMonnier et al. 



1999 



Massey et al. 



20061 ). A significant proportion of source 



D flux is missing from our observation when compared to single dish results (Figure [3]). 
The source we observe therefore appears to represent the warm core of a larger extended 
envelope found at the systemic velocity with lower average column density. This extended 



source is analogous to the s pherical wind descr i bed in previous mi 



obs ervations of VY CMa by 



and 



Tenenbaum et al. 



Ziurys et al 



fl2007h 



Ziurys et al 



fl2009h 



li meter wavelength 



Muller et al. 



fl2007h . 



(|2010[ ): however, only the first and last of these works observed SO2 
and did not identify it in the spherical wind. These previous works have found a relatively 
low expansion velocity of between 15 and 20 km s~^, which is high given the narrow velocity 
range of source in our channel maps (Figures H]- [7]). However, this discrepancy may be due 
to missing source D flux in our observations. The rotational temperatures derived from 
SMA and ARO data indicate the existence of a thermalized compact region with elevated 
temperatures with diameter 3 x 10^^ cm. In the sections below, we refer to this inner 
region as the core of source D. 

Our derived values of Trot and T^j^ distinguish between the lower temperatures of 
sources A, B, and C and a much hotter core of source D. In comparison to previous works, 



our temperatures 



Muller et al. 



'or so urce A, B, and C b racket the range of temperatures derived by 



fl2007h and 



Ziurys et al 



( 120091 ) (57K and 85K, respectively). This may be 
expected, as the preceding authors adopted the same best fit power-law temperature profile 
for both outflows, which are assumed to be at the same radius. As such, these previous 
results may represent an average of the temperatures of the outflows. 

Our inferred H2 densities fr om RADEX radiati ve transfer analysis are higher than the 



values found in previous studies. IZiurys et al 



(120091 1 adopted an isotropic II2 density profile 
based on an assumed mass-loss rate that gives a value of ~ 1 x 10^ cm^^ at 10^^ cm radius. 
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MuUer et al.l ( l2007l ) use a similar procedure to arrive at a lower value of 4.5 x 10^ cm ^ at 
the same location. A hig her density of 3 x 10^ cm~'^ is assumed for the inner wind region in 



Tenenbaum et al 



fl2007h . 



Part of this discrepancy between our values for H2 density and that of previous authors 
may be explained by the latter's assumption of isotropic mass flow, which does not account 
for density concentration in the spatially confined outflow regions. However, this reason 
alone may not be able to account for the more than two orders of magnitude difference. 
More significantly, our derived densities of between 5 x 10^ and 2 x 10^ cm^^ may be due to 
the presence of SO2 in regions of local density enhancement above the expected values from 
an isotropic model. We speculate that such high density regions are the result of shocks, 
and their existence around VY CMa can be inferred from the observations of OH masers, 



which overlap with sources A and B of SO2 emission ( iBowers et al. 



19831 ). The activation 

of the observed 1612 MHz maser line req uires H2 densities of between 1 x 10^ and 3 x 10^ 
cm~^ at 100 K (jPavlakis fc Kylafislll996l ). The upper end of this range is similar to our 



inferred value of minimum H2 density for sources B and C, while that of source A is much 
higher. However, OH masers may still be active in our source A despite its high densit y 



since its temperature of ~ 55K is cooler than that assumed in 



Pavlakis fc KvlafisI (119961 ). 



4.2. Sulfur chemistry in the CSE 



For the remainder of the Discussion, we address the implications of our observations 
for circumstellar chemistry by comparing SO2 distribution with t hose of the OH and 



SO molecules. 


SO^ ] 


Neiad & Millar 


1988 



Willacv fc Millar 



1997 



Cherchnefl 



20061 ): 



react ion ( IScalo fc Slavsky 



1980 : 



SO + OH 



SOo + H 



- 15 - 



The radial abundance of SO2 is therefore expected to reflect that of the two reactant 
molecules. A striking similarity between t he spectral profiles of SO2 and OH masers has 



already been noted by 



Ziurys et al. 



(120071 ). Maps of the 1612 MHz OH maser line show 
that it coincides with SO2 in sources A and B, but it is weak or undetectable in the outflow 
source C or the spherical wind sour ce D. The lack of O H maser detection in source C (or 



perhaps a very weak detection; see iBowers et al. 



19831 ) may be explained by anisotropic 



nature of maser radiation. Assuming that the velocity of outflows is oriented radially away 
from the star, relative velocities between different clumps of gas along a photon's line of 
travel are smallest when the path is parallel or antiparallel to the gas expansion velocity. 
The most efficient pumping of a masering state is then achieved when photons travel 
radially inward or outward from the star. T herefore, the strong est maser emissions are 



19761 ). This condition is nearly 



observed from sources along our line of sight (lElitzur et al.l 
met for sources A and B (15° to 22° from the line of sight). On the other hand source C, 
found at a line of sight angle of 54°, may also produce the 1612 MHz OH maser, but its 
signal is weak along the line of sight. 

The weak OH maser emission in the stellar velocity frame may be attributed to the high 
kinetic temperature of gas in this region. For a given volume density of gas, temperatures 
above a certain threshold tend to induce thermal equilibriu m in the emitting body, undoing 



the population inversions responsible for maser emissions (jPavlakis fc Kylafis 



19961). For 



temperatures of ~ 200K, the maximum allowable H2 density for the 1612 MHz OH maser is 
3 X 10^ cm~^, which is more than an order of magnitude lower than our inferred density for 
source D. The highly linear trend of Source D data in the rotational temperature diagram 
(Figure [9]) corroborates the prevalence of LTE conditions in source D. 

We therefore find that OH and SO2 distributions are generally similar and that their 
differences are reconcilable. 
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Finally, we ad dress the discrepancies between the SO2 and SO distributions of 



MuUer et al. 



( 120071 ). who have mapped the distribution of SO around VY CMa at a similar 



resolution to ours using a single rotational line of SO (J = 65 - 84; = 35K). SO was 
found in all four source regions described in this work, although the red-shifted SO emitter 
is not fragmented into two discrete sources as in the case ofSOo. Because of the strongly 



contrasting values of H2 density adopted in 



Muller et al. 



(I2OO7I ) and this work, direct 



comparisons between column densities are more instructive than comparisons between 
fractional abundances. For column density to act as a valid proxy for abundance, we must 
assume that the line-of-sight dimension of the corresponding SO and SO2 sources are similar 
and that the distribution of the each molecules within each source is homogeneous. While 
we cannot be certain that the line-of-sight thickness of the SO and SO2 sources are equal, 
their plane of sky dimensions are similar (~30% discrepancies). The combined column 
density of both SO outflows and the inner wind where they overlap along the line of sight 
was found to be ~ 10^^ cm~^. This value reflected the combined column density from both 
outflows and the spherical wind (corresponding to our sources A, B, and D). Given that the 
three regions contribute approximately equal amounts to the total SO column density, the 
value of Nso in each region is on the order of a few 10^^ cm~^. In contrast, we flnd column 
density of SO2 to be 1 to 2 x 10^^ cm~^ in each outflow source. The column density of SO2 
in the outflow sources A and C is therefore > 3 times greater than that of SO in the same 
regions. This statement likely applies as well to outflow source C, given its similar SO2 
column density compared to the other outflows and the lack of a discrete SO source at its 
location. Furthermore, even with missing flux, the column density of SO2 in the compact 
core of source D is within error as that of the outflow sources, implying that Nso2 > ^so 
in the core of the spherical wind < 3 x 10^^ cm from the star. 



Non-equilibrium C SE chemistry models ( IScalo fc Slavsky 



1980 



Willacv fc Millar 



Nejad fc Millar 



1988 



19971 ) make differing predictions about ratio of SO2 to SO abundance 
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{fso2/fso)- The IScalo fc SlavskyI ( Il980l ) (SS) model adopts a reaction rate for SO2 
formation that is fast compared to the rate of SO2 destruction via photodissociation. Under 
these conditions, SO is quickly converted into SO2, which becomes the primary S-bearing 
species. SO2 therefore reaches maximum abundance at a greater radius than SO. SS has 
predicted the radius of this transition where fso2/ fso > 1 to be around several times 10^^ 
cm. 

Our observations suggest that fso2/ fso > 1 for the outflow regions (> 4 x 10^® cm 
from star) and in the compact region within 3 x 10^^ cm of the star. If there does exist a 
transitional radius at which fso-,! fsn = 1, it must b e less than the latter radius, making 



Scalo fc SlavskyI (Il980l ) model. On the other hand, models 



our study consistent with the 
that predict fsn', / fsn < 1 at all radii are inconsistent with our data (jNejad fc Millar 



1988 



Willacv fc Millar 



19971 ). In addition, the SS model predicts a steady increase in the value 
of fso2l fso outward from the fso2l fso = 1 transition radius at ~ 2 x 10^^ cm. The value 
of fso2/fso increases by more than an order of magnitude by radius 10^^ cm. Observations 
of the outflow sources A, B, and C do not support this view, as fso2l fso in on the order of 
3 for all four sources. More precise comparisons between SO2 and SO abunda nce in each 



source is hindered by the lack of more detailed interpretations of SO maps (the 



Muller et al. 



( I2OO7I ) observation included only one SO line). Therefore, quantitative comparison between 
model radial distribution of the two species and our results is elusive and we regard our 
support of the SS model as only qualitative. 

Other uncertainties remain in our understanding of sulfur chemistry in the CSE of VY 
CMa. Our observations cannot be used to constrain whether the SO2 originates in the CSE, 



as assumed by the cited models , or in shocks in an "inner wind" wi thin a few ste 



the photosphere a s mod eled by ICherchnefB (120061 ) and observe d by 



Yamamura et al 



(I1999I ). Furthermore, as proposed earlier by 



Decin et al 



lar radii of 



Jackson fc Nguyen 



f|2010h a nd 



fll988h 
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and IWillacy fc Millarl ( 119971 ). shock events similar to those that occur in the inner wind 
may be the source of SO2 enrichment in the outflow lobes. Indeed, comparison of SO2 and 



maps of SO, PN, and SiS (IMuUer et al 



20071 : Figure [H this work) shows that SO2 exhibits 



a more fragmented distribution with discrete sources in the red-shifted outflow. Unlike 
these other molecules, the red-shifted SO2 outflow is partitioned into two regions of high 
local abundance (sources B and C), suggesting that local effects may participate in the 
fo rmation of SO^. H owever, CSE chemistry models that include the effect of shocks (similar 



to 



Cherchnefi 



20061) remain to be done. 



A further uncertainty involves the very high mass loss rate of VY CMa, which is, for 
exam ple, about two order s of magnitude faster than the asymptotic giant branch star IK 



Tau ( lOlofsson et al. 



19981 ). Varying the mass loss rate in CSE chemistry models does not 



signiflcantly affect the radial abundance proflles of chemicals species while it does result in 



Willacv fc Millaj fll997h . Therefore the 



globally larger envelopes for higher mass loss rates 
discussions of the f 302/! so proflle in this work are also relevant to stars with lower mass 
loss rates, while the actual measured radii of peak SO2 abundance around VY CMa are 
unique this object. 



4.3. Summary 

1. Four rotational lines of the SO2 molecule are mapped with ~ 1" resolution around 
VY CMa. SO2 is found in four discrete sources, three of which are fast (28 to 44 
km s~^) outflows far from the star and one is a slower spherical wind near the star. 
No symmetrical relationship among the faster outflows or visible and IR features are 
found. 

2. The three fast outflows are found at similar distances from the star and probably 
originated around 500 years ago. 
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3. Comparison between our SO2 maps and those of the 1612 MHz OH maser hne 
suggests that the two species are strongly correlated and that the OH maser detection 
may be limited by high temperature and density in the spherical wind. 



4. SO2 is more abundant than S O in all three outflow s ources, supporting the non- 



equilibrium chemistry model of 



Scalo fc Slavskv 



that predict fso2/fso < 1 for all radii (e.g., 



(119801). It is inconsistent with mode. 



Neiad fc Millar 



1988 



Willacy fc Millar 



19971 ). The distribution of SO2 in discrete clumps when compared to other molecules 



may point to the role of localized effects, such as shocks, in the enhancement of SO2 
abundance. 



We are grateful to Raymond Blundell, Thomas Dame and Patrick Thaddeus for 
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Table 1: Summary of observed SO2 lines. 



S02 


Frequency 




S 


transition 


(MHz) 


(K) 




42,2 — 3i,3 


235151.7 


19.0 


1.71 


103,7 - 102,8 


245563.4 


72.7 


5.44 


152,14 ~ 15l,l5 


248057.4 


119.3 


5.27 


16l,15 - 152,14 


236216.7 


130.7 


6.05 



E„ is the upper rotational state energy. 
"^S is the line strength. 



Table 2: Description of SO2 line emission sources* 



Source 


Offset 


Position 


Size 


Orientation 


{visr range in km s~^) 


(") 


Angle (°) 


(") 




A: (-18, +10) 


0.66 + 0.10 


33° 


2.5 X 2.4 


126° 


B: (+26, +66) 


0.83 + 0.10 


316° 


2.9 X 1.3 


41° 


C: (+26, +66) 


2.45 + 0.09 


272° 


2.9 X 1.8 


127° 


D: (+12, +26) 


~ 


N/A 


3.2 X 1.9 


40° 



* Elliptical gaussian fits of the identified SO2 emission sources. Offsets are measured in arcseconds from the 
center of the continuum emission (RA=07''22™58.336^ Dec=-25°46'03.063"). Source size represents the 
major and minor axis of the ellipse, and its orientation angle is that of the major axis. Angles are given as 
rotations from north through east. Fits were performed on maps of the 235151.7 MHz line. 
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PN J=5-4 




4 ^ 2 ^ ^ ^2 ^ ^ 
E— W position offset (orcseconds) 



Fig. 1. — Left: PN J=5-4 and SiS J=13-12 spectra showing the peak emission at 19.49 km 
s~^ (from the hne rest frequencies). We adopt this velocity as that of the stellar frame. Right : 
Integrated intensity maps of PN (top) and SiS (bottom) emission, which show maxima close 
to the peak of the continuum emission. 
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1 1 1 1 1 1 1 

'■ 16(1,15)-15(2,14) 
^236216.6 MHz 

[ j 


1 1 1 1 1 1 1 1 . 
1 130 K ■ 


:i 1 1 1 1 1 1 

'■ 15(2,14)-15(1,15) 
^248057.4 MHz 


119 K ■ 

-| 1 1 1 1 1 1 H 


H It ^1 i 1 1 1— 

'■ 10(3,7)-10(2,8) 
^245563.4 MHz 


73 K ; 

\ ' ^ 


J 1 1 1 1 1 1 

: 4(2,2)-3(l,3) 
-235151.7 MHz 


1 19 K ; 

1 1 " 



-100 100 

V. Velocity (km s"') 



236216 
I I I I I I I I I 




248057 







T 




245563 
^1 I I I I I I I I 




a> 
u 



OO - 



u 
a> 
o 



-5 



235151 
I I I I I I I I I 




5 0-5 
RA Offset 



Fig. 2. — Left: SO2 spectra averaged in a region of 3"x3" around the continuum peak. 
Right: Integrated intensity maps corresponding SO2 lines over velocity intervals —20 and 
10 km s~^ are shown in blue contours. The red contours show the integrated intensity over 
velocity intervals 20 to 60 km s~^. The starting value and interval of contour levels are set 
to 35 mJy beam~^ km s"^. The continuum emission is shown in grey scale (repeated in all 
panels). The location of the three identified sources A, B and C is shown in the bottom two 
maps. 
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-60 -40 -20 20 40 60 80 100 

VIsr velocity (km/s) 



Fig. 3. — Left: SO2 spectra from our (black) and 10 meter single dish observations by the 
Arizona Radio Observatory's Submillimeter Telescope (red) for the two lowest temperature 
transitions (235151.7 and 245563.4 MHz). The SMA data was smoothed at the appropriate 
scale to approximate the field of view of the single dish telescope ( 30"). 
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Fig. 4.— Channel maps of SO2 4(2,2)-3(l,3) emission at 235151.7 MHz toward VY CMa. 
In each panel, emission is integrated over a 10 km s^^ wide velocity range centered on the 
velocity indicated on the top right corner. Contour levels are same as in Fig. O 
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Fig. 5.— 



Same as in Figglfor SO2 10(3,7)-10(2,8) emission at 245563.4 MHz. 
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Fig. 6.— 



Same as in Figglfor SO2 15(2,14)-15(1,15) emission at 248057.4 MHz. 
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h. , I , . . . 




5 -5 5 -5 5 -5 

RA Offset 



Fig. 7.— Same as in Figglfor SO2 15(1,15)-15(2,14) emission at 236216.6 MHz. 
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Fig. 8. — Left: A position velocity cut along east-west at Declination offset of 0" in the SO2 
line emission at 245563.4 MHz, at position angle 105°. The contours correspond to 10% of 
the peak emission. Sources A, B, C and D are clearly separated. 



- 32 - 




140 



Fig. 9. — Rotational temperature diagram of the four SO2 emission sources. Red, blue, and 
black curves and points indicate sources A, B, and C respectively. Orange represents the 
spherical wind (source D). Variable "L" on the vertical axis represents the left-hand side 
of Equati on [H The dashed o range line represents the spherical wind based on ARO single 



dish data 



Tenenbaum et al. 



( l2010h . The errors plotted are set to three sigma. Horizontal 



coordinates of the plotted points have been staggered for clarity. 
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Fig. 10. — Fits between our measured line ratios and predictions with Radex radiative 
transfer. Solid, dashed, and dotted contours enclose regions of test p-values greater 
than 5%, 30%, and 95%. Owing to lower S/N and weak dependence of line ratios to these 
quantities, the derived temperature and density of source D are much less well-constrained 
than the other three. 
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Table 3: Summary of LTE and RADEX modeling results 


Source 


Trot 


SO2 Column Density 


Tkin 


H2 Density 


SO2 Abundance ( /502) 




(K) 


(cm~^) 


(K) 


(cm~^) 




A 


611?? 


9.2ljj X 10^^ cm-2 


50 


> 2x10^ 


7 X 10"^° 


B 


noli 


1.9l[i| X 10'' 


~ 75 


> 3x10^ 


2 X 10-^ 


C 


691?^ 


X IQis 


- 100 


> 3x10^ 


8 X 10^9 


D 


2401^^ 


7.6 ± 1.0 X 10^5* 


- 220 


> 4x10^ 


5 X 10-1° 



* Values for source D suffer from missing extended flux. Using single dish data, average T^ot = 97 K and 
N = 1.8 X 10^^ cm~^ assuming that the emitter is an uniform 4" by 4" region. 



